Normal modes of relativistic systems 
in post Newtonian approximation 



Y. Sobouti 1,2 ' 3 and V. Rezania 



1. Institute for Advanced Studies in Basic Sciences, P. 0. Box 45195-159, Zanjan, Iran. 
OO ■ 2. Physics Department, Shiraz University, Shiraz, Iran. 

as ; 

OS . 3. Center for Theoretical Physics and Mathematics, AEOI, Tehran, Iran. 
$_i ' e-mail: sobouti@sultan.iasbs.ac.ir and rezania@sultan.iasbs.ac.ir 




> 



^1" 

o 

OO 
OS 



Publication in main journal 



a,; 



Proofs to be sent to V. Rezania 



To appear under section 5, stellar cluster and associations 



o 

H 



CO 



> 




1 



Abstract. We use the post Newtonian (pn) order of Liouville's equation (pnl) to study the 
normal modes of oscillation of a relativistic system. In addition to classical modes, we are able to 
isolate a new class of oscillations that arise from perturbations of the space-time metric. In the 
first pn order; a) their frequency is an order q smaller than the classical frequencies, where q is 
a pn expansion parameter; b) they are not damped, for there is no gravitational wave radiation 
in this order; c) they are not coupled with the classical modes in q order; d) in a spherically 
symmetric system, they are designated by a pair of angular momentum eigennumbers, (j,m), 
of a pair of phase space angular momentum operators (J 2 , J z ). Hydrodynamical behavior of 
these new modes is also investigated; a) they do not disturb the equilibrium of the classical 
fluid; b) they generate macroscopic toroidal motions that in classical case would be neutral; c) 
they give rise to an oscillatory component of the metric tensor that otherwise is zero for 
the unperturbed system. The conventional classical model are, of course, perturbed to order q. 
These, however, have not been discussed in this paper. 
Key wordrrelativistic systems: oscillations- normal modes. 

1 Introduction 

In recent years some authors have debated the existence of a new class of quasi modes in rela- 
tivistic systems believed to have been originated from the perturbations of the space-time metric. 
Kokkotas and Schutz (1986) considered a toy model of a finite string (to mimic a fluid) coupled 
to a semi-infinite one (to substitute the dynamical space-time). Such a system accommodates 
a family of damped oscillations. Baumgarte and Schmidt (1993) obtained the same results by 
considering a slightly more realistic toy model. Kojima (1988) verified that strongly damped 
(w-) modes do indeed exist in realistic stellar models. Subsequently, Kokkotas and Schutz (1992) 
used a WKB method to study w mode spectra for several polytropic models. Different inves- 
tigators have proposed different mathematical and numerical schemes to isolate these modes. 
Discrepancies exist. Leins, Nollert and Soffel (1993) employed two different approaches for the 
modes of Schwarzschild black hole. They used Leaver's (1985) continued fraction approach and 
a Wronskian technique (Nollert and Schmidt 1992). They confirmed the findings of Kokkotas 
and Schutz. They, however, found other modes with smaller oscillation frequencies and stronger 
dampings than those that had been found by the WKB scheme. Anderson, Kokkotas and Schutz 
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(1995) used a numerical approach to clarify the reason for the discrepancies. They came to agree- 
ment with the findings of Leins et. al. (1993), and showed that the WKB method employed by 
Kokkotas and Schutz failed to distinguish between these modes. Recently, Lindblom, Mendell 
and Ipser (1997) have used a formulation of the relativistic pulsation equation in terms of two 
scalar potentials, one for the fluid and the other for the perturbations in the gravitational field. 
They have found a set of even-parity modes for a relativistic star, including the w modes. 

In this paper, we approach the problem through relativistic Liouville's equation rather than 
relativistic hydrodynamics. The reason for doing so, is to avoid thermodynamic concepts incor- 
porated into hydrodynamics. Liouville's equation is a purely dynamical theory and free from 
such complexities. One, of course, pays the price by having to dwell in the six dimensional 
phase space an elaborate mathematical task, but not obtrusive. In section 2 we give the pn 
order of the linearized Liouville equation that governs the evolution of small perturbations from 
an equilibrium distribution function. In section 3 we extract the evolution equation for the non 
classical modes and explore some of its properties. In section 4 we explore the 0(3) symmetry of 
the modes and classify them on basis of this symmetries. In section 5 we study hydrodynamics 
of the pn modes. In section 6 we on a variational approach to the solution of pn modes and 
give a numerical analysis for the modes of polytropes. Lengthy mathematical manipulations are 
given in appendices A and B. 



2 Formulation of the problem 

Liouville's equation in the post-Newtonian approximation pnl for the one particle distribution 
of a gas of collisionless particles maybe written as follows 

H^ + £ d + (7 £ P >(x,u,t)=0, (1) 

where (x, u) are phase space coordinates, q is a small post-Newtonian expansion parameter, the 
ratio of Schwarzchild radius to a typical spatial dimension of the system, say. The classical and 
post-Newtonian operators, C cl and C pn , respectively, are 

c " = -<»'s* + £?57>- < 2 *> 

„ P n ... 9 , ns d9 a a 89 a89 8® A,drii drij, dn, 8 ,„ , 
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The imaginary number i is included for later convenience. The potentials 0(x,t), 0(x, i) and 
r/(x, i), solutions of Einstein's equations in pn approximation, are 

9(x , t) = /£^V , (Xlt) = 4/i!^|fV (So,*) 

where dT = d 3 xd 3 u. See Rezania and Sobouti (1998, hereafter paper I) for details. In an 
equilibrium state, F(x, u) is time-independent. If, further, it is isotropic in u, macroscopic 
velocities along with the vector potential 77 vanish. It is also shown in paper I that the following 
generalizations of the classical energy and classical angular momentum are integrals of pnl: 

e = e d + qe pn = 1^ -e + q{26 2 -Q), (4a) 
k = EijkX-i u k exp{q6) lf(l + q9), for spherically symmetric 6{r) and 0(r). 

(4b) 

Equilibrium distribution functions in pn approximation can be constructed as appropriate func- 
tions of these integrals. In paper I the pn models of polytrope were studied in this spirit. 

Here we are interested in the time evolution of small deviations from a static solution. Let 
F — > F(e) + 6F(pc,u,t), \ 5F |<C F V(x, u, t). Accordingly, the potentials split into large and 
small components, 0(r) + 89(x,t), 0(r) + <50(x, t) and 5n(x,t) where r = |x|. Both, the large 
and small components, can be read out from Eqs. (3). Substituting this splitting in Eq. (1) 
and keeping terms linear in SF gives 

+ C cl + qC pn )5F + {5C cl + q5C pn )F(e) = 0, (5) 

where C cl and C pn are now calculated from Eqs. (2) with 0(r), 0(r) and n = 0, and SC cl and 
5C pn with 50(x,t), <50(x,t) and 5n(x,t). For the latters, Eqs. (2) give 

6C * F = -i^^ = F e C d 59, (6a) 
dx l du l y ' 

5C pn F = F e [£ cl (5@ - 4956) +iu 2 ^-- iu 1 ^}, (6b) 

where F e = dF(e)/de. The classical limit of Eq. (5) is 

Hi- + C cl )SF + F e C cl 59 = 0, S6 = -G f ^''"'^ rfT'. (7) 
at J x — x' 



Equation (7) is the linearized Liouville-Poisson equation. It was studied briefly by Antonov 
(1962). He separated 6F and Eq. (7) into even and odd components in u, extracted the 
equation governing the odd component and showed that 5F 0( id satisfies an eigenvalue equation. 
Sobouti (1984, 1985, 1986, 1989a, b) elaborated on this eigenvalue problem, studied some of 
its symmetries and approaches to its solution. Sobouti and Samimi (1989), and Samimi and 
Sobouti (1995) showed that Eq. (7) has an 0(3) symmetry and its oscillation modes can be 
classified by the eigennumber j and m of certain phase space angular momentum operator. 

The effect of the pn corrections in Eq. (5) on the known solutions of the classical Eq. (7) 
can be analyzed by the usual perturbation techniques. Whatever the procedure, the first order 
corrections on the known oscillation modes will be small and will not change their nature. We 
will not pursue such issues here. The main interest of the this paper is to study a new class 
of solutions of Eq. (5) that originate solely from the pn terms and have no precedence in 
classical theories. It is not difficult to anticipated the existence of such modes. Functions of 
the classical energy and angular momentum satisfy Eq. (7) with no time dependence. They 
may be considered as the eigenfunctions of that equation belonging to the infinitely degenerate 
eigenvalue to = 0. Now the pn corrections of Eq. (5) set up an eigenvalue equation in this 
degenerate subspace and give rise to a new class of modes that have relativistic origin and are 
not coupled to the classical solutions, in the first pn order. 

3 The post-Newtonian modes of oscillations 

Theorem: Any function G(e d , \ d ) of classical energy and angular momentum, e cl and l cl , respec- 
tively, for which 5/o(x) = / Gd 3 u = 0, is a static solution of the linearized Liouville-Poisson, Eq. 
(7). For, C d G = and <50(x) = f 6p(x!)\x. - x^dV = 0, QED. 

One class of such functions are those with odd parity in u or, equivalently, in 1, i. e., 
G(e cl ,l d ) = -G(e d , -l d ). We will see in Eqs. (17) and (18) below that there exists a class of 
odd parity functions that not only fulfill the requirements of the theorem but also give 5@ = 0. 
With these preliminaries, we come back to solutions of Eq. (5). 

Let 5F(x,u,t) = 5F(e d ,lf)e- iu}t , where 5F(e d , if) is a static solution of Eq. (7). By the 
theorem above, Eqs. (5) and (6) reduce to 

q£ pn 5F = oj(5F + qF e u i 5ri i ). (8) 
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Furthermore the term containing 5r) is an order q smaller than 8F and can be dropped. Thus 

qC pn 8F = uj5F. (9) 

where 

O 7 

C pn = —{[(u 2 -4e)9' + Q']x i -e'(ic-u)u i }^-, '"" = -?-, (9a) 
r ait* ar 

see Eq. (A. 1) for details. We owe the reader a demonstration of 8Q = 0. In appendix A we 
show that oj is real. There follows that 5F is complex; for, C pn is purely imaginary. Thus, 
5F = 5F R + iSFj . Equation (9) can be decomposed accordingly: 

qC pn 8F R =ioj8F n (10a) 
qC pn 8F I = —iuj8F R . (10b) 

Upon eliminating the imaginary component one obtains 

q 2 C pn2 8F R =co 2 8F R . (11) 

The same equation is satisfied by 5F r The latter will, however, be obtained from Eq. (10a) once 
8F R is calculated from Eq. (11). In appendix A we show that C pn2 is Hermitian, in spite of the 
fact that C pn is not. Furthermore, C pn2 is even in both x and u. Therefore, the eigensolutions 
SF R can be chosen either odd or even in x and u. Even parities, however, are not acceptable. 
For they give rise to non vanishing changes in the macroscopic density, 8p, and the potential 89 
and 8Q, which violates the validity of Eq. (8). The x and w-parities of 8F n however, will be 
opposite to those of 8F R ; for they are connected by Eqs. (10a, b), where C pn is odd in x and u. 

4 0(3) symmetry of D> n 

In appendix A we show that L pn and C pn2 commute with the phase space angular momentum 
operators Jj = —i£ijk{x 3 + u^-^). Thus, they can have simultaneous eigensolutions with 
J 2 and J z . Eigensolutions of the latter pair of operators, Aj m ; j and m integers, are worked 
out in appendix B. They are complex polynomials of the components of the classical angular 
momentum vector. Their x and u parity is that of j. 

An expression of the form / (e c/ , l d2 ) A jm is already an eigensolution of J 2 and J m belonging 
to eigennumbers j and m. See appendix B, Eqs. (B. 2-5). Solutions of Eq. (11) for 8Fr will, 
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therefore, reduce to substituting an expression of this sort with an odd j value in Eq. (11) and 
solving it for f(e cl , l cl2 ). In section 6 we will take a varitional approach to such solutions, and 

as variational trial functions will consider functions of e cl only. For the moment let 

N 

SF R = f jm = f(e)Re A jm = [ ^ c n (-e) n ]ReA jm , c n = consts, (12) 

n=j+l 

be a solution of Eq. (11) belonging to the eigenvalue u>. Combining this with its corresponding 
imaginary component from Eq. (10a) we obtain 

8F jm (x, u, t) = (1 + Id") fjne-™. (13) 

id 

This expression satisfies Eq. (9), if fj m satisfies Eq. (11). At this stage let us note an important 
property of Liouville's equation. If a pair (oj, SF) is an eigensolution of Liouville's equation, 
(— uj,5F*) is another eigensolution. This can be verified by taking the complex conjugate of 
Eq. (9). These solutions, being complex quantities, cannot serve as physically meaningful 
distribution functions. Their real or imaginary part, including time dependence, however, can. 
With no lost of generality we will adopt the real part. Thus, 

Re 5Fj m (x, u, t) = f(e)Re A jrn cos tot + i-C pn (f(e)Re A jrn ) smut. (14) 

This, however, is the sum of two eigenfunctions belonging to oj and — oj and is not a simple 
eigensolutions of Eq. (9). 

5 Hydrodynamics of pn modes 

In this section we calculate the density fluctuations, macroscopic velocities, and the perturba- 
tions in the space-time metric generated by a pn mode. 

For j an odd integer, /j m (x, u) is odd while C pn fj m is even in x and u. The macroscopic 
velocities are obtained by multiplying Eq. (14) by u and integrating over the u-space. Only the 
odd component of 5Fj m contributes to this bulk motion, 



pv 



j f(e)Re Aj m ud 3 u cos ut. (15) 



In appendix B we show that pv is a toroidal spherical harmonic vector field. In spherical polar 
coordinates it has the following form 

—1 d dY 

p(v r , v#, Vtp) = r J G(v es )(0, Re £^Q^ Y 3rn{#,<p), Re -^(&,tp)) cosut, (16a) 



where 



Jo 

and v es = y/29 is the escape velocity from the potential 0(r). The macroscopic density, generated 
by even component of Eq. (14), is 




6p(x,t) = i- [ C pn (f(e)Re A jm )d\smujt 




The second integral is obtained by an integration by parts. The vanishing of it comes about 
because of the fact that the radial vector x is orthogonal to the toroidal vector pv. One also 
notes that V • (pv) = 0. It can further be verified taht, the continuity equation is satisfied at 
both classical and pn level. To complete the proof of validity of Eq. (9) one should also show 
that SO = 0. From the definition Eq. (3c) for 50 and Eq. (14) for 5F one has 



The first two terms are zero because 5p and 56 are zero. The third term vanishes because the 
integral over u' has the same form as in 5p except for the additional scalar factor u /2 in the 
integrand. Like 5p it can be reduce to the inner product of the radial vector x and a toroidal 
vector that are orthogonal to each other. QED. 

The toroidal motion described here slides one spherical shell of the fluid over the other 
without perturbing the density, the Newtonian gravitational field and, therefore, the hydrostatic 
equilibrium of the classical fluid. In doing so, it does not affect and is not affected by the 
conventional classical modes of the fluid at this first pn order. 

Nonetheless, the pn modes are associated with space time perturbations. From Eq. (8c) of 
paper I and Eq. (3b) of this paper, goi component of the metric tensor is 




(18) 




(19) 



In spherical polar coordinates, one has 



rjr = 0, 



(20a) 
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1 d 

= -ajRe - — z-z-Yjm^, ip) cosui, (20b) 
sm v op 

dY 

7]^, = ajRe ^ m (■&, p) cos ut, (20c) 



where 



16vr 



(r/Ry yj (R) + (2j + 1)H J r R r'- j - 1 y j (r')dr' forr<R 
(R/ry +1 yj (R) forr>R 



J 2j + 1 

y ,(r) = r- J ' _1 / ' r /2j ' +2 G(e(r / ))dr / , (20e) 
Jo 

G(0(r)) = f(e)u j+3 du 
Jo 

= 2 j / 2+1 T(j/2 + 2)r(n + l)0(r) n+J '/ 2 + 2 /r(n + j/2 + 3), (20f) 

where R is the radius of the system and T(n) is the gamma function. The remaining components 
of the metric tensor remain unperturbed. 

6 Variational solutions of pn modes 

In Eq. (12) we suppress the subscript jm, for brevity. It will be remembered that / will be 
chosen odd in both x and u and will be a function of the phase coordinates through the classical 
energy and angular momentum integrals. We multiply Eq. (11) by /* and integrate over the 
volume of the phase space available to the system 



lo 2 



J fdF = q 2 J fC pn2 fdT = q 2 J (C pn ^fy(£ pn f)dT. (21a) 



^ J fdF = J \C pn f\ 2 dT + J 6' 2 [u 2 - 3^-^]/ 2 dr. (21b) 



The adjoint operator C pn ^ is given in Eq. (A. 5). Substituting the latter in Eq. (21a) and simple 
manipulations gives 

^ [ fdF = I \£ pn f\ 2 dF + f 6' 2 \u 2 - 3 (X ' U) ^ 

It will be remembered that oj is real and oj 2 is positive. Therefore, in spite of one negative term, 
the right hand side of Eq. (21b) has to be positive. That is, £P n2 is positive in addition to being 
Hermitian. To solve Eq. (21b) by variational techniques we choose the set of trial functions 
{(— e) n Aj m , j = odd, n > j, e = u 2 /2 — 6}. Any member of this set is an eigenfunction of 
J 2 and J z with eigenvalues j and m, respectively. Next we choose a real linear combination of 
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these trial functions as an ansatz for / 

N 

fjm = f{e)Re A jm = [ c n (-e) n ]ReA jm , c n = consts. (22) 

n=j+l 

We substitute Eq. (22) in (21b), carry out the integrations, and minimize oo 2 with respect to 
variations of c n 's. We arrive at the following matrix equation to solve for oj 2 and c„'s: 

uj 2 

WC = -^SC, (23) 



where C is the column vector of the variational parameters, {c n }, and elements of the S and W 
matrices are 

S pq = J (-ef+^Re A jm | 2 dr, (24a) 

= J (CP n (-e)PRe k jm )*{CP n {-e)iRe k jm ) + j 9' 2 [u 2 - 3^-^-]{-e)^\Re A jm \ 2 dT. 

(24b) 



w pq 



Eigen w's are the roots of the characteristic equation 

> 2 _ 

,2 



W - ^S\ = 0. (25) 
Q 



For each oo, Eq. (23) can then be solved for the eigenvector C. This completes the Rayleigh-Ritz 
variational formalism of solving Eq. (11). In what follows we present some numerical values for 
polytropes. 



pn Modes of polytropes belonging to (j, to) = (1,0) 

From Eq. (B. 9), A\ q = l z = rncos7, where 7 is the angle between x and u. Substituting this 
in Eqs. (24) and integrating over directions of x and u vectors and over < u < y/26 gives 

S - f 1 6 p+q+2 - 5 x 4 dx, (26a) 
Jo 

W pq = 47rGp c {(16a pq -b pq ) / 9' 2 9 p+q+3 - 5 x 4 dx 

Jo 

+(l-8a pq ) ( ®'0'6 p+q+2 *x A dx 
Jo 
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+a pq f 1 e' 2 9 p+q+1 - 5 x 4 dx}, (26b) 
Jo 

_ pqjp + q + 2.5) _ 4p + 4g + 11 23 ... f26c 

Pq ~ (p + q)(p + q-l)' Pq ~ P + q + 3.5' P ' Q ~ ' ' ' ( 



Polytropic potentials 9 and 6 were obtained from integrations of Lane Emden equation and Eqs. 
(28) of paper I, respectively. Eventually, the matrix elements of Eqs. (26), the characteristic Eq. 
(25) and the eigenvalue Eq. (23) were numerically solved in succession. Tables 1-4 show some 
sample calculations for polytropes 2, 3, 4, and 4.9. Eigenvalues are displayed in lines marked 
by an asterisks. The column following an eigenvalue is the corresponding eigenvector, i. e. the 
values of ci, C2, • ■ of Eq. (22). To demonstrate the accuracy of the procedure, calculations 
with six and seven variational parameter are given for comparison. The first three eigenvalues 
can be trusted up to four to two figures. Convergence improve as the polytropic index, i. e. the 
central condensation, increases. Eigenvalues are in units of AitGp c q 2 and increase as the mode 
order increases. 
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Appendix A: Some properties of C pn and C pn2 

The Hilbert space: let 7i be the space of complex square integrable functions of phase coordinates 
(x, u) that vanish at the phase space boundary of the system, 

H : {#(x, u), J g*gdT = finite, g{boundary) = 0}. 

Integrations in Ti are over the volume of phase space available to the system. In particular, the 
boundedness of the system in velocity space sets the upper limit of u at the escape velocity, 
y/29, where #(x) is the potential at x. Thus, g(x, \/26{x)) = 0. 

The operator C pn is defined in TL. For a spherically symmetric, isotropic and static system, 
Eq. (2b) gives 

C pn = Z!{[( n 2 _ 40 j 9 > + Q>y _ 0/( x . u)^}^-. (A.l) 

We have set 77 and all time derivatives equal to zero, used d/dxi = (xi/r)d/dr and denoted d/dr 
by "/ " . 

Theorem 1. Eigenvalues of C pn are real: This is in spite of the fact that £ pn is not Hermitian. 
Proof: Let 

C pn f = <jf, f € H. (A.2) 
Noting that C pn is purely imaginary, the complex conjugate of this equation is 

C pn f* = -oj*f*. (A.3) 

Multiplying (A. 2) by /* and (A. 3) by /, adding the two equations and integrations over the 
phase space gives 

= J (f*c pn f + fC pn f*)dT = J £ pn {f*f)dT 

= (-2*) J e'^ffdr, 

= (u- W*) J ffdT. (A.4) 

The x and u parities of /*/ are even and those of x • u are odd. Therefore, the integral vanishes 
and u) = uj* = real. QED. 

Hermitian adjoint of C pn : For f,g£Tt, by an integration by parts one finds: 

J g*C pn fdT = J[(C pn + 2z^e')g]*fdr, 
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thus 

X • u 



C pn\ = jyn + 2i — 9'. (A.5) 

r 

Theorem 2. C pn2 is Hermitian in Ji. Proof: For g G Ji 
J g*C pn2 gdT = J (C pn ^g)*C pn gdr 



[(C pn + 2i 9'g\*C pn gfT 



r 



-I 

= J \C pn g\ 2 dT + j 6' 2 {u 2 - 3 (X r2 U) ]g*gdT = real, (A.6) 

QED. Some integrations by parts are carried out in calculating the term containing x • u. 
0(3) symmetry of C pn and C pn2 : For a spherically symmetric system, C pn of Eq. (A. 1) 
depends on the angle between x and u and their magnitudes. Simultaneous rotations of the x 
and u coordinates about the same axis by the same angle should leave C pn form invariant. The 
generator of such simultaneous infinitesimal rotations on the function space H is 

J t = j} = -is ijk (x^ + A). (A.7) 

With 9 and G of Eq. (A. 1) being spherically symmetric, it is easy to show that C pn commutes 
with Jj 

[£ pn , Ji] = 0. (A.8) 

The operators Ji, however, obey the angular momentum algebra. Therefor, the three operators 
C pn , J 2 and J z commute pairwise. The same is true of C pn2 , J 2 and J z . In appendix B we 
elaborate on the eigenfunctions of J 2 and J z to prepare the path for simultaneous eigensolutions 
of the latter trio. 
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Appendix B: Eigensolutions of J 2 and J z 

As pointed out earlier, Jj's of Eq. (A. 7) have the angular momentum algebra, 

[Ji, Jj] = iEijkJk- (B.l) 
The simultaneous eigensolutions of J 2 and J z , A jm (x, u), obey the followings 

J 2 A jm = j(j + l)Ajm, j = 0, (B.2) 



«/ z A im = mA jm , -j <m< j. (B.3) 
The ladder operators, J± = J x ± iJ y , raise and lower m values: 



J±A jm = yj (j Tm)(j±m + l)A jm± i. (B.4) 

In particular 

J±A j: ±j = 0. (B.4a) 

The effect of J, L on classical energy integral, e = u 2 /2 — 0(r), and the classical angular momentum 
integral, k = e follows 

J,e = Jil 2 = Jif(e,l 2 ) = 0, (B.5a) 

Jilj = isijkh- (B.5b) 

It is simple to verify Eqs. (B. 5) by direct substitution. 
Theorem 1: 

A,. : , 4 !•'(/, : (B.6) 

Proo/: 

J*4 = jl 3 ± l (JJ±) = ±jl{, by (B.5b), (B.7a) 

J 2 /^ + = (j_ J + + J 2 2 + J 2 )/^ = j(j + 6 y (5.4a) and (5.7a), (B.7b) 

j 2 zi = ( J+J _ + j2 _ = j ( i + (b.7c) 

QED. Combining Eqs. (B. 6), (B. 4) and (B. 5) one obtains 

A jm = af(e, l 2 )J{~ m P_ = bf(e, l 2 )Jr m l{, (B.8) 
14 



where /(e, I 2 ) is an arbitrary function of its arguments, and a and b are normalization constants. 
Examples: 

Ai = k, (B.9a) 

Ai±i = /±, (B.9b) 

A 20 = 2/+/_-/2 = I(3/ 2 2 -/ 2 ), (B.9c) 

A 2 ±i=Z±Z*, (B.9d) 

A 2 ±2 = l±. (B.9e) 

Theorem 2: The vector field V- 5 "™ = / Aj m udft is a toroidal vector field belonging to the spherical 
harmonic numbers (j,m), where integration is over the directions of u. 

Preliminaries: Let ip) and (a, (3), and 7 denote the polar angles of x, of u and the angles 
between (x, u), respectively. Also choose magnitudes of x and u to be unity, as only integrations 
over the direction angles are of concern. One has cos 7 = cos?? cos a + sini? sin a cos(c^ — (3) 

u r = cos 7, (B.lOa) 

u$ = — sin?? cos a + cosi? sin a cos (99 — /?), (B.lOb) 

Utp = — since sin(</? — (B.lOc) 

Z+ = z(sin# cos a e llf - cos?? sin a e l/3 ). (B.lOd) 

Proof: By induction, we show that a) V" is a toroidal field and b) if V jm is a toroidal field, so 
is V'™" 1 . 

a) Direct integrations over a and [3 gives 

V™ = /" l j + u r dn = 0, dfi = sin a da d(3, (B.lla) 
^ = / liu#dn = -^°- Yjj (i>,<p), (B.llb) 
V£ = J l{ Uip dn = Ay..^^). q^d. (B.llc) 

b) Suppose V- 7 " 1 is a toroidal vector field and calculate V J m_1 = J(J_Aj m )udfl, where J± = 
L± + K±, L± = ±e ±iLp (§ 5 ±icotg°d-^) 1 K± = ±e ±i(i (-^±icotga-^). Again direct integrations 
gives 

Vj m " 1 = L_y/ m = 0, if V r jm = 0, (B.12a) 

15 



m-1 1 8 , r / „ n . „ T ,1 m 13 



^ m "=— rjr-Yjm-^,?), if Vj™ = - — —Yj m {0, <p), (B.12b) 

^'^yj^p), t/ ^ m = ^ m (^). (B.12c) 



QED. 
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Table 1: pn mods of polytrope n=2, belonging to (j, m) = (1,0). Eigenvalues are in units 
4nGp c q 2 . c n 's are linear variational parameters. The first three eigenvalues are reliable 
up to three figures. Higher order eigenvalues are less accurate. A number a x 10 ±fe is 
written a a ± b. 

Table 2: Same as Table 1. n — 3 and (j, m) = (1, 0). 
Table 3: Same as Table 1. n = 4 and (j, m) = (1,0). 
Table 4: Same as Table 1. n — 4.9 and (j, m) = (1,0). 
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Table 1. 



cu 2 


.1825+01 


.4973+01 


.6448+01 


.1216+02 


.3425+02 


.1686+03 




Cl 


.3113+02 


-.8912+02 


.1663+03 


.1344+03 


.7545+01 


-.1399+04 




c 2 


.3908+02 


.1045+04 


-.3234+04 


-.9746+03 


-.2392+04 


.8484+04 




c 3 


-.1420+03 


-.6649+04 


.1801+05 


.4514+04 


.7952+04 


-.9647+04 




c 4 


.5803+03 


.1804+05 


-.4351+05 


-.7014+04 


-.2607+03 


-.2251+05 




c 5 


-.9110+03 


-.2210+05 


.4724+05 


.8324+03 


-.1811+05 


.5188+05 




c 6 


.5252+03 


.1020+05 


-.1874+05 


.2882+04 


.1317+05 


-.2717+05 




oo 2 


.1823+01 


.4865+01 


.5895+01 


.9113+01 


.1465+02 


.4228+02 


.3226+03 


Cl 


.3028+02 


-.7086+02 


.1529+03 


-.3129+02 


.1561+03 


-.4624+02 


.2042+04 


c 2 


.4812+02 


.6908+03 


-.2810+04 


.1313+04 


-.1513+04 


-.2762+04 


-.1461+05 


c 3 


-.1305+03 


-.3993+04 


.1702+05 


-.5686+04 


.6685+04 


.1077+05 


.2271+05 


c 4 


.2576+03 


.8181+04 


-.4788+05 


.3425+04 


-.3673+04 


.1875+04 


.4154+05 


c 5 


.1303+03 


-.3086+04 


.6823+05 


.2433+05 


-.2910+05 


-.4718+05 


-.1496+06 


c 6 


-.7534+03 


-.7924+04 


-.4771+05 


-.4855+05 


.5132+05 


.5873+05 


.1425+06 


c 7 


.5475+03 


.6707+04 


.1302+05 


.2568+05 


-.2386+05 


-.2120+05 


-.4423+05 



pn x pn 2 pn 3 pn 4 pn 5 pn 6 pn 7 
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Table 2. 



o 

u 2 


.1534+01 


.4836+01 


.9473+01 


.1938+02 


.4083+02 


.1128+03 




Cl 


.9752+02 


-.6975+02 


.2464+03 


-.2246+03 


-.9102+03 


.3169+04 




c 2 


.3284+02 


-.8725+03 


-.1121+04 


-.2590+04 


.1713+05 


-.2631+05 




c 3 


.2096+03 


.3859+04 


.5591+04 


.1444+05 


-.1023+06 


.6390+05 




c 4 


-.5354+03 


-.5728+04 


-.1216+05 


-.9903+04 


.2599+06 


-.3406+05 




c 5 


.3941+03 


.2528+04 


.5215+04 


-.2221+05 


-.2933+06 


-.4814+05 




c 6 


.1803+01 


.1125+04 


.3307+04 


.2153+05 


.1208+06 


.4268+05 




oo 2 


.1533+01 


.4688+01 


.7993+01 


.9068+01 


.1124+02 


.1909+02 


.1093+03 


Cl 


.9318+02 


-.1440+03 


-.1202+03 


-.1069+04 


-.5706+03 


-.5482+02 


.3703+04 


c 2 


.1121+03 


.6997+03 


.5482+04 


.1856+05 


.7685+04 


-.5626+04 


-.3381+05 


c 3 


-.2118+03 


-.4506+04 


-.2955+05 


-.1063+06 


-.4112+05 


.3078+05 


.1007+06 


c 4 


.2709+03 


.9777+04 


.5298+05 


.2726+06 


.7791+05 


-.4371+05 


-.1109+06 


c 5 


.1206+03 


-.9309+03 


-.6283+04 


-.3375+06 


-.1278+05 


-.7049+03 


.1239+05 


c 6 


-.7005+03 


-.1574+05 


-.7154+05 


.1894+06 


-.9027+05 


.3228+05 


.4581+05 


c 7 


.5309+03 


.1200+05 


.5087+05 


-.3511+05 


.5945+05 


-.1218+05 


-.1722+05 



pn x pn 2 pn 3 pn 4 pn 5 pn 6 pn 7 
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Table 3. 



cu 2 


.7569+00 


.2822+01 


.5661+01 


.8814+01 


.1519+02 


.6952+02 




Cl 


.6291+03 


-.1067+04 


.2143+04 


-.1949+04 


-.6870+04 


.1400+05 




c 2 


-.9217+02 


.1770+04 


-.1693+05 


.1131+05 


.8373+05 


-.2337+06 




c 3 


.4162+03 


.2808+04 


.5682+05 


-.3654+04 


-.3195+06 


.1293+07 




c 4 


-.3883+04 


.5860+04 


-.1184+06 


-.2807+05 


.4791+06 


-.3112+07 




c 5 


.6427+04 


-.2303+05 


.1257+06 


-.4668+04 


-.2545+06 


.3371+07 




c 6 


-.3089+04 


.1612+05 


-.4514+05 


.3416+05 


.1251+05 


-.1344+07 




oo 2 


.7569+00 


.2813+01 


.5021+01 


.8747+01 


.1272+02 


.3322+02 


.7683+02 


Cl 


.5590+03 


-.8716+03 


.2653+03 


-.2421+04 


.1881+04 


.1412+05 


.3376+05 


c 2 


.1189+04 


-.2018+04 


.1406+05 


.1926+05 


-.7436+04 


-.2356+06 


-.5191+06 


c 3 


-.6377+04 


.2349+05 


-.1057+06 


-.4732+05 


-.5363+05 


.1298+07 


.2528+07 


c 4 


.9376+04 


-.3509+05 


.2059+06 


.6165+05 


.2228+06 


-.3112+07 


-.4750+07 


c 5 


.5449+03 


-.4645+04 


-.2977+05 


-.4272+05 


-.7106+05 


.3356+07 


.2298+07 


c 6 


-.1192+05 


.4364+05 


-.2533+06 


-.3854+05 


-.4046+06 


-.1333+07 


.2455+07 


c 7 


.7228+04 


-.2275+05 


.1775+06 


.5845+05 


.3227+06 


-.1382+03 


-.2085+07 



pn x pn 2 pn 3 pn 4 pn 5 pn 6 pn 7 
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Table 4. 



cu 2 


.4481+00 


.1827+01 


.4078+01 


.6515+01 


.1170+02 


.1391+03 




Cl 


-.2888+02 


.1663+03 


-.2794+03 


.1593+03 


.1405+03 


.1081+05 




c 2 


-.2440+03 


-.7593+04 


.2050+05 


-.2099+05 


.2665+05 


-.2129+06 




c 3 


.4933+05 


-.2772+04 


-.1400+06 


.1883+06 


-.3467+06 


.1344+07 




c 4 


-.1722+06 


.1443+06 


.2902+06 


-.5138+06 


.1372+07 


-.3583+07 




c 5 


.2124+06 


-.2675+06 


-.2194+06 


.4871+06 


-.2092+07 


.4207+07 




c 6 


-.8916+05 


.1394+06 


.5712+05 


-.1179+06 


.1073+07 


-.1790+07 




oo 2 


.4380+00 


.1805+01 


.4006+01 


.6190+01 


.7980+01 


.1439+02 


.8964+02 


Cl 


-.1701+02 


.1379+03 


-.3341+03 


.3427+03 


-.3020+03 


.7695+03 


.8642+04 


c 2 


-.6649+03 


-.6322+04 


.2326+05 


-.3097+05 


.2196+05 


-.1111+05 


-.1534+06 


c 3 


.5135+05 


-.1143+05 


-.1601+06 


.2940+06 


-.2552+06 


.1349+06 


.8174+06 


c 4 


-.1667+06 


.1599+06 


.3264+06 


-.9227+06 


.1022+07 


-.9712+06 


-.1565+07 


c 5 


.1694+06 


-.2551+06 


-.1784+06 


.1132+07 


-.1574+07 


.3018+07 


.4968+06 


c 6 


-.1770+05 


.8656+05 


-.9582+05 


-.4879+06 


.7432+06 


-.3959+07 


.1421+07 


c 7 


-.3646+05 


.3341+05 


.9586+05 


.2938+05 


.8318+05 


.1819+07 


-.1047+07 



pn x pn 2 pn 3 pn 4 pn 5 pn 6 pn 7 
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